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We consider the dynamics of lattices which have constrained constitutive units flexible 
in only their mutual orientations. A continuum description is derived through which it 
is shown that the models have zero shear velocity, free-particle like internal rotational 
modes and volume decreasing linearly with temperature. The relevance of models to a 
range of problems is pointed out. 



INTRODUCTION 

The dynamical properties of solids with strong anharmonicity are often usefully studied by specifying 
the constraints on the atoms rather than the forces between them—This approach has yielded interesting 
results in phenomena ranging from the properties of interfaces!!! to criteria for glass-formationH'tl to 
rigidity percolationH in solids with " floppy Modes" . 

We consider here the problem of lattices in which some (strongly anharmonic) bonds can be deformed 
at a much lower energy than the others. As noted below the problem is of wide interest but our attention 
to it is drawn by the contraction in volume of a class of solids linearly as temperature is increased, m*& 
range of temperature which in some cases extends from ~ 10K to the highest measured temperaturdj'tl 
Conventional perturbative calculations with anharmonic potentials are not very helpful in this case. 
These solids have a complicated lattice structure but always contain polyhedra which share some or all 
corners with the neighbors. Any distortion of the polyhedra is expensive, the vibrations changing the 
angles between the polyhedra however have low energies ~ loq. A simple model can be constructed to 
show the phenomena of thermal contraction in such cases. A two-dimensional version is shown in fig. 
la, which consist of a solid with three squares per unit cell, two of which are connected only at two 
corners (these can be replaced by sticks of fixed length) while the third is connected at all four corners!! 
All bond-lengths as well as the internal right-angles of the squares are constrained to be fixed. The only 
independent variable is the angle Oi between the squares. 

The area occupied by the solid is the sum of the area occupied by the solid squares and the area 
between them. The latter are bounded by perimeters whose length are fixed. Under this condition the 
area is least when the perimeter has the highest symmetry, i. e. in the configuration shown in fig. la. 
Any variation in {Oi] must then necessarily lead to a reduction in volume. For thCpjiiiit cell shown, the 
area decreases as Of. We expect (0 2 ) ~ T for T > cj , giving the observed behaviorp'El. 

The above is the correct dimensional argument for the thermal contraction, but it should be clear that 
's in different cells are coupled, and coupled non-locally to the strain. The coupling of translational 
and rotational degrees of freedom leads in the continuum limit to a field-theoretical problem with some 
remarkable features. Besides the thermal contraction, we have found that the shear modes have zero 
stiffness in such solids and that they have rotational modes with dispersion characterestic of free particles: 
u> oc fc 2 . These features are also present in three-dimensional versions of the model and are expected also 
to hold in more general models than studied here. 

A convenient set of variables for the problem are the position of the four corners of each square i, 
{\Pi, n = R-i + r i,n} * with Ri the coordinate of the center of the square i and n = ±1, ±2 label each 
corner (see Fig. lb). Fixing the lengths and the internal right-angles of the squares, the relatives 
positions r n j can be expressed as functions of the angle 6i alone, r^„ = b(cos(0i + (n — l)ir/2)x + 
sin(8i + (n — l)n/2)y), where & is a fixed length. The original 8 degrees of freedom per unit cell are then 
reduced to 3 {uj, Oi}, where U; — — R° = (ui X x + Ui V y) is the displacement of the center of mass of 
the squares around the equilibrium position at T = 0. Two more constraints per unit cell come from 
fixing the lengths of each connecting stick (or square); 

|*i,i - *i+ x ,-i\ = |*i,a - *i+»,- 2 | - 26 (1) 
the labels ii x, and iiy represent the neighbors of the cell i. 



In any distortion about the symmetric position, the two-diagonals in each square change angle with 
respect to the connecting sticks . The potential energy may be expressed in terms of the four such angles 
per unit cell, 



V = K ^2 C 1 - C0S Mi,n)) 



(2) 



i,n=±l,±2 

fii^n may be written in terms of the vectors ^j, n . For example 

(# M - *i _i) • (*i+x -1 - 
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In the continuum limit a9 Q r <C 1, a9 Q u < 1 (a= 46, a = (a;, y)), the potential energy simplifies to 



(3) 



y = 



l 



n • (24 + ^u) + r 2 • (2y + <9 y u)) + i(r x ■ fign + r 2 • flgr 2 ) 



Note that terms quadratic in 9 Q u do not occur in Eq. (g). This Equation can be written as: 



2 - cos $ (2 + Vu) -sin^ (V x u) - 2b 2 {V 6)' 



dS. 



(4) 



(5) 



The rotational invariance of Eqs. (g) and (||) will be evident only after the constraints of Eq. ([!]) are 
implemented. 

The kinetic energy is T — \m^ li n (*$>i.n) 2 which in the continuum limit becomes 



T=^ I (M(u) 2 +Amb 2 {9) 2 )dS 



(6) 



To obtain the continuum version of the constraints of Eq.(l), consider all the four sticks connecting 
any given square to its neighbors. The four equations can be written as 



|A Q | 2 + |B Q | 2 ±2A a -B Q = 46 2 
or \K a \ 2 + \B a \ 2 =46 2 
A Q • B Q = 



(7) 



where, 



A Q = a(& + <9 Q u)-2r-— d 2 r 

a 2 

B a = yd 2 u-a<9 Q r, 



(8) 



These four equations determine the four components of the strain tensor d a u. We have found that 
these constraints are satisfied to 0(9 2 ) by two sets of relations, as may be verified by substitution 
Set I: 



Set II: 



a(x + d x u) = 4ri + -a 2 d 2 ri, 
3 

a(y + d v u) = 4r 2 + ja 2 d 2 r 2 . 



a(x + d x u) = 2(n + r'O + yf x + a 2 (^ 2 n + ^d 2 x r[) 
a(y + d y u) = 2(r 2 + r 2 ) - xf y + a 2 {\d 2 y v 2 + ^d 2 r' 2 ) 



(9) 



(10) 



where r' n = b(cos(8 + (n — l)n/2)x — sin(d + (n — l)7r/2)y)and f x (9) and f y (0) are arbitrary real 
functions of order 9 2 or higher. These relations relates the deformation of the unit cell with 9. For Set 
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I the first term on the right side rigidily rotates each unit cell by 8, a uniform 9 gives a rigid rotation of 
the entire solid. The second term on the right hand allows a smooth spatial variation of 9 (see Fig. 2). 
The deformation for the set II are shown in Fig. lb. In these optical-like modes the connecting stick 
rotates by —9 and the cell contracts ~ 9 2 in each direction. 

Now we determine the equation of motion for each set of solutions. Eqs. ([)[) can be combined in the 
form 



V • u = 2(cos6» - 1) - 3& 2 V.V(cos<9) = J {9), 
V x u = 2sin(9 + 3fr 2 V.V(sin(9) = J x (0). 



(11) 



3q(9) and Ji(0) are thus the sources determining the irrotational and solenoidal parts of u. The other 
two equations obtained from Eqs.Q are that the shear strain (d x u y + d v u x ) — 36 2 (<9 2 — d 2 )(sin9), and 
(d x u x — dyUy) — 3b 2 (d 2 — d 2 )(cos9). For results at long wavelengths these can be equated to zero. The 
model thus has zero shear-velocity. 

With Eqs. (H), it is straight forward to show that 



V 



-4Kb 2 



\V9\ z dS. 



(12) 



The absence of terms proportional to 9 in 
traced to the local invariance satisfied by dO) 



(|T2|) has important consequences. Their absence may be 



V x u 



V • u 



36 2 |V6>r = 1 



(13) 



Corrections to Eq. ^ are of 0(6 4 (V6») 4 ). 

Consider next the equation of motion for 9(r, t). To do this u in the kinetic energy T must be 
expressed in terms of the sources J (0), Jj.(0). Noting that |J (6»)| 2 + |Ji(6»)| 2 = 49 2 + 12b 2 9V 2 9 + 0(9 3 ), 
the equation of motion for 9 to linear order associated with set I is found to satisfy the energy-wave 
vector relation 



(m - 6M)lu 2 - AKq 2 + AM = 



(14) 



The last term comes from (u) 2 . Eq. ( |l4|) yields the interesting result that at long- wave-length, free- 
particle like modes exist for which u —jjq 2 . The physical reason for the existence of such rotational 
mode lies in the local invariance ([jj]) and the fact that rotation of the entire solid (q = 0) involves 
transporting macroscopic mass. At long- wavelengths the effective mass is ~ M/q 2 . The characterestic 
frequency below which the dispersion of the modes is proportional to q 2 , determined by looking at 
corrections of 0(q 4 )) in Eq.@, is O(KfM) 1 / 2 . 
For Set II the potential energy is given by 



V = AK I [1 - cos 29 (1+b 2 |V0| ) + sin 9 (f x - f y )] dS 

and the equations for the sources of the fields for V.u and V x u are, respectively 

J o (0) = 2 cos 9 -2-36 2 (cos 9 (V6») 2 + sin6» V 2 6»), 
■MO) =f x -fy + b 2 (cos9 \7 2 9 - sin9 (V<9) 2 ). 



(15) 



The contribution of the source-fields to the kinetic energy is now ~ 9 A and does not change the 
equation of motion for 9(r, t) which give the expected optical-like dispersion relation 



8K b 2 q 2 
M 2 '' 



(16) 



To linear order f a do not enter in the dispersion. 

Since the problem has only one degree of freedom per unit cell, the free-particle like and the optical- 
like modes are not linearly independent. This may be seen through Fig. 2. The free-particle like modes 
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can be locally thought of as composed of a rotation by 9 + 59 and a q — it /a optical mode of amplitude 
69/2. The sharing of the spectral weight between collective modes in non- linear problem is known in 
other contexts. 

Now we come to the question of thermal compression with which we started. The equilibrium volume 
of the lattice at a temperature T may be determined from the thermal expectation value 

(V.u) = (Jo(0)>. (17) 

From Eqs. ( |TT| ) and ([l5|), we note that both sets of modes contribute to the thermal contraction at all 
temperatures which to the order of the theory here is proportional to the expectation value — (9 2 (T)). For 
T larger than 0(ojq) = 0(K /M) 1 / 2 the thermal contraction may then be determined by the equipartition 
theorem to be proportional to T. The thermal contraction is entropic akin to that of rubber. 

The problem can be extended to 3 dimensions replacing the squares by octahedra and adding an 
extra stick per cell which connects the corners of neighboring octahedra in the z direction. The model 
has a cubic symmetry. For a fixed size of the octahedra the 6 degrees of freedom per unit cell (3 
translation and 3 rotations) are reduced to 3 by fixing the length of the sticks (analogues of Eq.(l)). 
the coordinates of the 6 vertices of each octahedron can again be expressed as *&i,n = Ri + r- ln , with 

n,i = ~ri,-i = (6, 0,0)C W , r i)2 = -rj,_ 2 = (0, b, 0)0^ and r[ 3 = -v\ » = (0,0, b)C^ where 
C^e^> = CI CgC^ represent the three rotations by the Euler angles <f>, 9 and ijxA. 

The 3 constraints in the continuum limit can be expressed in the same form as Eq. (Q), where now 

a 2 

A a = a(e a + d a u) - 2r„ Q — ^d 2 r Ua 
a 2 

B a = —d 2 a u-ad a r na . 

A a and B Q are three component vectors, a — (x,y,z), e a is the unit vector in the a direction and 
n a = (1, 2, 3), indicate the corresponding corners of the octahedra. Once again the shear velocities are 
zero. Also we again get two sets of solutions. 
Set I: 

3 

a(e a + d a u) = 4r„ Q + -a 2 9 2 r„ Q (18) 



Note that r„ Q are functions of the Euler angles, then again eqs. (|18|) give the deformation of each unit 
cell as a function of the Euler angles. 

The potential energy can be found analogous to the 2D case (Eq. Q2h and results 



V 3D = -b 2 8K 



Tr C7 B i, ■ V 2 Q 



dV (19) 



with (s7 2 C ( j,0^ — d 2 (Cj,8ti>)i.j- The absence of a mass term is again due to an equations analogous 

to Eqs. @. We again obtain 3 modes with uj ~ q 2 for reasons similar to the 2-dimensional case. 
For Set II 

a(e Q + d a u) = 2(r fe . + r' fc .) + a 2 (±d 2 x r na + \d 2 x v' n J, (20) 

where is an arbitrary vector satisfying \r' n \ = 2b. From here it can be derived the three optical modes. 

We have also considered the simple the model in which the sticks connecting the rigid polyhedra 
(squares in 2-dimension) have finite elasticity, i.e. the constraints of Eqs. (0) are not imposed. The model 
then has of course the usual longitudinal and transverse modes at long wavelengths. The results obtained 
above cannot be obtained in any limit of the parameteres of such harmonic model as a pertubative 
expansion about them. The constraints lead to qualitativaly new class of properties. These properties 
are preserved if wegeneralize the sticks in Fig. (2) to a polymer, i.e. with a chain of fixed length sticks 
deformable in the mutual orientation at the connecting points. 

The model investigated is especially interesting in some unfamiliar contexts. It has direct application 
to sintered solidstJ, where rigid particles are connected to each other through bonds with small energies 
for angular variations although the model needs to be generalised to take into account the random 
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placements of the bonds. The predictions for zero or very small shear velocity and of long wavelength 
rotational modes is expected to hold generally and should be tested in that case. 

Some models for granular solidstHl are related to the model studied here again with generalisation to 
the random placing of the bonds. The prediction of zero shear velocity obtaiped here is particularly 
interesting in connection with granular materials because they flow like a liquid. L3 Of-course to examine 
flow the model must also be generalised to allow for " bond-breaking" . The model studied here is also a 
particular example of general models discussed by S. Alexander^ with potential applications to unusual 
excitation modes in glasses and foams. 
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FIGURE CAPTIONS 

Figure 1 The unit cell of the model of rigid squares hinged at corners, (a) At the minimun of 
potential, (b) with a given distortion exhibiting compression in the optical mode. 
Figure 2 exhibits the distortions in the rotational mode. 
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